function [S,HL,EdX,SDdX] = wrapper_counterfactual(pars,M_full,W_full,Nu,verbose)     
    
output                                  = evaluate_cf(pars,M_full,W_full,Nu,0);
X_data                                  = output.OutputFinal.X_monthly;
[T,Nd,Ns]                               = size(X_data);
X_data                                  = reshape(X_data(:),[T Nd*Ns]);
dX                                      = (X_data(end,:) - X_data(6,:))';
S                                       = round(100*pars.S,1);
HL                                      = round(-log(1-0.50)/pars.theta,1);
EdX                                     = round(100*mean(dX),1);
SDdX                                    = round(100*std(dX),1);
if verbose
    fprintf('Shock size  (pp)                   % 3.1f \n',S);
    fprintf('HL          (months)               % 3.1f \n',HL);
    fprintf('E(dX)       (pp)                   % 3.1f \n',EdX);
    fprintf('sd(dX)      (pp)                   % 3.1f \n',SDdX);
    fprintf('\n');
end

end